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• The optimization model minimizes annual biomass-ethanol production costs. 

. The model optimizes strategic and tactical decisions simultaneously. 

• Miscanthus-ethanol production costs in a base-case were $0.72 L -1 of ethanol. 

• A 50% reduction in biomass yield would increase unit production costs by 11%. 

• Higher demand would reduce unit production costs by achieving economies of scale. 
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To ensure effective biomass feedstock provision for large-scale biofuel production, an integrated biomass 
supply chain optimization model was developed to minimize annual biomass-ethanol production costs 
by optimizing both strategic and tactical planning decisions simultaneously. The mixed integer linear 
programming model optimizes the activities range from biomass harvesting, packing, in-field transpor¬ 
tation, stacking, transportation, preprocessing, and storage, to ethanol production and distribution. The 
numbers, locations, and capacities of facilities as well as biomass and ethanol distribution patterns are 
key strategic decisions; while biomass production, delivery, and operating schedules and inventory mon¬ 
itoring are key tactical decisions. The model was implemented to study Miscanthus-ethanol supply chain 
in Illinois. The base case results showed unit Miscanthus-ethanol production costs were $0.72 L _1 of eth¬ 
anol. Biorefinery related costs accounts for 62% of the total costs, followed by biomass procurement costs. 
Sensitivity analysis showed that a 50% reduction in biomass yield would increase unit production costs 
by 11%. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Cellulosic based biofuels are considered to be sustainable 
renewable transportation fuel. To meet a target of 16 billion gal¬ 
lons of cellulosic ethanol production in 2022 (EPA, 2010), more 
than 200 million Mg of biomass will be required annually consid¬ 
ering a cellulosic ethanol conversion rate of -70-90 gallons per Mg 
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of dry cellulosic biomass (Wyman, 2007; Humbird et al„ 2011). 
Few commercial cellulosic ethanol facilities, however, exist due 
to a lack of a cost-effective technology and reliable feedstock sup¬ 
ply. Bioprocessing facility related costs and biomass procurement 
costs comprise the major biofuel production costs (Humbird 
et al„ 2011; Lin et al„ 2013). Biomass is produced in a distributed 
manner within a limited harvesting window each year; but it 
should be processed in a centralized facility throughout the year 
to gain the economies of scale. Therefore, effective and efficient 
provision of biomass is a key challenge for large-scale biofuel 
production. 
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Supply chain management has been used to facilitate effective 
biomass provision and, generally, three levels of decision-making 
are considered: strategic, tactical, and operational decisions. We 
seek to focus on the interplay between strategic and tactical 
decision-making. Strategically, biomass resource evaluation and 
selection of facility location and capacity are important long-term 
decisions. Biomass availability and associated production costs 
have been evaluated for feedstocks including Miscanthus (Khanna 
et al., 2008), switchgrass (Perrin et al., 2008), and corn stover 
(Kadam and McMillan, 2003). Furthermore, considering the spatial 
variances of biomass availability and production costs, several 
models have been presented to optimize facility locations and 
capacities (Kim et al., 2011; Lin et al., 2013). These strategic plan¬ 
ning models provide decision support based on annual biomass 
delivery estimations, without much consideration for tactical plan¬ 
ning details. 

Tactically, biomass can only be harvested within a limited win¬ 
dow due to its standing dry matter loss. Biomass production, deliv¬ 
ery, and operating schedules and inventory monitoring are key 
tactical decisions. Biomass procurement consists of multiple unit 
operations including harvesting, packing, in-field transportation, 
and handling. Several biomass supply chain models have been 
developed to manage biomass production and delivery activities 
for different feedstocks, including Miscanthus (Shastri et al., 
2010), switchgrass (Kumar and Sokhansanj, 2007; Zhu et al., 
2011; Sharma et al., 2013), wood residues (Svanberg et al., 2013), 
corn stover (Sokhansanj et al., 2006; Leboreiro and Hilaly, 2011), 
and cotton stalks (Tatsiopoulos and Tolis, 2003). These models, 
however, are based on the given strategic decisions such as deter¬ 
mined facility locations and capacities. 

Strategic decisions of biomass supply chain will impact subse¬ 
quent tactical decisions. Without the support of biomass delivery, 
processing facilities cannot achieve their designed operating 
capacity. Few multi-scale supply chain optimization models have 
been developed to solve processing facility locations and capacities 
and biomass delivery schedules simultaneously (Eksioglu et al., 
2009; Zhang et al., 2013). However, these studies did not consider 
much on farm management for biomass production. 

Determining how to harvest, store, and deliver biomass to sup¬ 
port the processing activities is vital to optimizing biomass supply 
chains. Therefore, it is important to coordinate and optimize link¬ 
ages between biomass production, logistics, processing simulta¬ 
neously. This study aims to develop an integrated biomass 
supply chain optimization model to minimize annual biomass- 
biofuel production costs. The model simultaneously optimizes 
both strategic decisions, such as facility locations, capacities, and 
resource allocations, and tactical decisions, such as biomass pro¬ 
duction, delivery and operating schedules as well as inventory 
monitoring. 


2. Methods 

2.1. Integrated model overview 

Large-scale biofuel production faces significant challenges 
including logistics, delivery schedules, and inconsistent feedstock 
formats (Hess et al„ 2007). To facilitate effective biomass produc¬ 
tion, a five-stage biomass supply chain was proposed to include 
farms, centralized storage and preprocessing (CSP) facilities, biore¬ 
fineries, ethanol blending stations, and ethanol consumption. 
Farmers were expected to harvest, bale, transport, and stack bio¬ 
mass at their farm gates for delivery. CSP facilities were expected 
to store, handle, and preprocess biomass, in order to provide con¬ 
sistent feedstock format for conversion and reduce logistics chal¬ 
lenges (Hess et al., 2007; Eranki et al„ 2011; Lin et al„ 2013). 


After production at biorefineries, ethanol will be transported to 
blending stations to mix with gasoline. The mixed fuel will be 
transported to gas stations for end consumption. For this study, 
we assume current blending infrastructure and have been well 
developed and no new blending facilities will be built. Therefore, 
the objective of the integrated model is to optimize a five-stage 
biomass-biofuel supply chain with the focus on the development 
of biomass production, storage, preprocessing, and conversion 
infrastructure and the selection of existing ethanol blending 
infrastructure. 

By the nature of biomass production, spatial and temporal con¬ 
straints are important for biomass supply chain optimization. Spa¬ 
tially, number of farms and cropland size, and biomass yield vary 
by county. The transportation distance between facilities varies 
based on their spatial distribution. These spatial features are highly 
related to strategic decisions such as the selection of facility loca¬ 
tions and capacities. Temporally, the harvestable biomass yield 
and probability of worldng day vary with time. Using these tempo¬ 
ral constraints tactical decisions, such as when to harvest, store, 
and deliver biomass, can be considered. Without coordinating 
these activities, the loss of biomass in operations could offset the 
gains of the high harvestable yield. 

To address spatial and temporal issues, the integrated optimiza¬ 
tion model is composed of four modules: a biomass farm manage¬ 
ment module that optimizes biomass harvesting, packing, loading, 
and infield transportation activities; a logistics planning module 
that optimizes biomass delivery schedules and transportation 
fleet; a location-allocation module that optimizes facility locations 
and capacities, transportation modes, and biomass distribution 
patterns; and an ethanol distribution module that optimizes etha¬ 
nol distribution patterns (Fig. 1). The model coordinates tactical 
decisions, such as production and delivery schedules, with strate¬ 
gic decisions such as facility locations and capacities, to provide 
decision support on both levels accordingly. This coordination 
was implemented by linking delivery schedules among biomass 
production activities at farms, biomass utilization activities at 
CSP and biorefinery facilities, and ethanol demand at consumption 
areas (Fig. 1). 

The integrated optimization model is a mixed integer linear 
programming (MILP) model that was developed on the GAMS plat¬ 
form (GAMS, 2012) and was solved using the CPLEX® solver. A list 
of set names, decision variables, and parameters used in the model 
is provided in “Nomenclature” in Appendix. The model was devel¬ 
oped through the integration of key equations from the previously 
developed tactical (Shastri et al., 2010) and strategic models (Lin 
et al., 2013) to coordinate and optimize the entire supply chain 
as a single problem rather than considering strategic or tactical 
decisions individually. The integrated model reduces some unnec¬ 
essary assumptions used in the optimization, considers the inter¬ 
actions of decisions between two levels, and provides a high 
resolution and quality of decision support accordingly (Table 1). 
In addition to being a MILP model that optimizes strategic spatial 
decisions such as farm and facility locations, the integrated model 
is also a discrete time optimization model that optimizes tactical 
temporal decisions such as biomass harvesting and delivery sched¬ 
ules. The objective of the integrated model is to minimize annual 
biomass-ethanol production costs (Z), which consists of five parts: 
biomass procurement costs ( Z f ), transportation costs (Z v ), CSP facil¬ 
ity related costs (Z p ), biorefinery related costs (Z r ), and ethanol dis¬ 
tribution costs (Z d ) (Eq. (1)). 

Z = Z f + Z v + Z p +Z r + Z d (1) 

2.1.1. Farm management module 

Biomass production on farms consists of four major unit opera¬ 
tions: harvesting, packing, in-field transportation, and handling. 
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Fig. 1. The components of the integrated optimization model and their data flow. The rounded rectangles represent source data (cropland area and farm number (USDA, 
2010), biomass yield Gain et al., 2010), land opportunity costs (Khanna et al., 2008), probability of working day (Shastri et al., 2012), equipment and vehicle data (Shastri et al., 
2010), road networks (US Census Bureau, 2011), and preprocessing and storage data (Hess et al., 2007)). The rectangles represent the models used in the system 
(MISCANMOD (Jain et al., 2010), NREL model (Humbird et al., 2011), and Shortest Distance Estimation Model and Facility Cost Estimation Model (Lin et al., 2013)). Farm 
management, logistics planning, facility allocation, and ethanol distribution modules are the models developed in this study. The clouds represent output data. 


Table 1 

The improvement of the integrated model for decision support can be illustrated by considering the variation in its resolution, key assumptions and inputs to the models, and key 
decision variables and model outputs relative to that of the previously developed strategic and tactical models. Notice that several key assumptions and inputs used previously by 
the strategic and tactical models are now notable outputs from the integrated model. Please note, additional inputs and assumptions are necessary to operate the integrated 
model, but listed here are the key differences in the three models. 


Strategic model (Lin et al., 2013) 


el ( r al., 2010) 


County level 

-Biomass availability and production 
costs at each county 
-Constant monthly delivery schedule 
-An estimation of biomass loss rate 
in the system 


-Preselected farms, processing facility 

locations and capacities 

-Preselected annual biomass flow pattern 


Key Decision Variables and -Optimal facility locations and 

Model Outputs capacities 

-Optimal annual biomass flow 
patterns 


-Optimal selection of farm equipment 
and transportation vehicle 
-Optimal biomass harvesting and delivery 
schedules 

-Biomass availability and production 
-Biomass loss rate throughout the system 


-Biomass availability and production costs at 
each county 

-Optimal selection of farms and processing 
facility locations and capacities 
-Optimal selection of farm equipment and 
transportation vehicles 
-Optimal biomass harvesting and delivery 

-Biomass loss rate throughout the system 


Biomass procurement costs are estimated from the equipment cap¬ 
ital related costs, operating costs, and opportunity costs of crop¬ 
land allocated to energy crops (Eq. (2)). Given the nature of crop 
growth and standing loss, biomass needs to be harvested in a lim¬ 
ited time period. If more equipment can be operated during the 


high biomass yield period, more biomass could be harvested, 
which could reduce the unit biomass operating costs; however, 
this would increase equipment capital related costs. Different 
types of equipment have their unique operating capacity and 
performance that could be suited for different sizes of farms. 
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Therefore, farm selection (0} J ), equipment selection (F'/ m> ), and 
harvesting and other farm operating schedules (p/' m ‘' ) are the 
key decisions for farm management, where i represents biomass 
supply county, j represents farm number, t represents time period, 
and m f represents farming equipment (e.g. harvesting and baling 
machines). Related to these decision variables, annualized equip¬ 
ment cost {//'), unit operating costs farm area (a iJ ), and 

county-level cropland opportunity costs (z' 0 ) are the key input 
parameters. 

% = EEE F / ,mf >< V +EEE Df* x 

+ EE a ' J x 0/ x T„ (2) 

The harvesting and other farm operating schedules (P'/' m,t ) are 
on a daily basis. It is assumed that these schedules remain constant 
during the time period, which typically consists of several days. 
The size of farm affects the biomass production costs, where large 
farms could purchase and implement large equipment gaining 
high productivity; however the BioFeed model cannot select the 
appropriate farm for production (Shastri et al., 2010). Farm selec¬ 
tion is a key strategic decision added in this integrated model, to 
help biofuel producers procure biomass efficiently (Eq. (3)). The 
harvesting schedule is designed to monitor the cumulative har¬ 
vested area of each farm by that time period (ff 1J t ). If a farm is 
selected for biomass production (the binary variable 0/ will be 
set equal to 1), all its farmland area (ac' J ) needs to be harvested 
completely by the end of the harvesting season (t = t h ); whereas 
if a farm is not selected, its cropland will not be harvested at all. 

H i# * = a iJ x 0/ (3) 

The equations regarding farm equipment selection [F /' m ’) and 
biomass production schedules are based on the previously 

developed BioFeed model (Shastri et al., 2010). The model assumes 
that all harvested biomass needs to be packed, transported, and 
stacked to be available for delivery at farm gate within the same 
time period. Therefore, given the amount of biomass required for 
operating during each time period, the type and number of equip¬ 
ment required for each unit operation can be optimized. Given the 
optimized harvesting schedules and the temporal changes of bio¬ 
mass yield, the amount of daily biomass delivered to farm gate 
at any time period (B| Jt ) can be estimated. Detailed constraints 
regarding farm activities can be found in Shastri et al. (2010). 

The harvested biomass inventory at each farm at the end of 
each time period {I/' 1 ) is monitored based on the accumulation of 
the differences between the daily biomass delivered at farm gate 
and its daily delivery schedules to all the CSP facilities ( D ,Jk,t ) while 
considering daily biomass losses at farms (a)/), where A t represents 
the length of that time interval in days (Eq. (4)). Since biomass 
inventory cannot be negative at any time, the amount of biomass 
delivered from each farm is limited by biomass availability at the 
farm gate, which links farm management with logistics planning. 
The model assumes all the biomass at each farm needs to be deliv¬ 
ered to the CSP facilities by the end of each year (t= f) (Eq. (5)). 

fU = - E^'*) x A t j x (1 - ay x A t ) (4) 

t/^ = 0 (5) 

2.1.2. Logistics planning module 

Biomass transportation costs (Z v ) consist of fleet equipment 
capital related costs [Z Vc ), operating costs (Z„J, and handling costs 
(Z Vh ) (Eqs. (6)-(9)). The model assumes road transport using 


truck-trailers as the only option. The model assumes the biofuel 
producers own the transportation fleet and can manage logistics 
to realize just-in-time delivery for processing. Key decisions in¬ 
clude delivery schedules (T LJkm " x and T k,l,m " ,t ) and fleet vehicle 
selections (F™"). Related to these decision variables, annualized 
vehicle costs (/""), unit operating costs (Py") and handling costs 
(/ 4 t and p hi ) are key input parameters. 


z v =z Vc +z Vo +z Vh 

(6) 

z Vc = Y/7 x K” 

(7) 


=(EEEEE^ W x ^ +EEEE T ' I, ' , ' T “" t x 

x/C 

(8) 

^ = EEEE° yAt x n + EEE 0 ^ x /4, (9) 

Weight and volume constraints are the two major issues for 
biomass transportation. Different types of vehicles have different 
weight (V™") and volume capacities (1/™^,). Based on weight and 
volume constraints of each type of vehicle and biomass density 
data (<7 iJ ), the number of trips required for each type of vehicle 
during each time period (T' J ’ k ’ m "’ t ) is estimated (Eq. (10) and 
(11)). Similar equations will be applied to estimate the number 
of trips needed for each type of vehicle between CSP facilities 
and biorefineries (T klm, ' x ). The shortest travelling distances be¬ 
tween facilities are calculated based on the existing road networks 
(Lin et al„ 2013). Given the typical travelling speed, the shortest 
travelling time between facilities for each type of vehicle ( d ,j,k,mv 
and is estimated. Given the planned delivery schedules 

and travelling time as well as daily work time (r/), the daily number 
of vehicles needed during each time period (N™ vX ) is determined 
(Eq. (12)). To ensure there exist sufficient vehicles for delivery at 
any time period, the total fleet size should not be smaller than 
the required number of vehicles at any time period (Eq. (13)). 

D im < x T 1J,k ’ m "’ t (10) 

c«) 

K"* x n > EEE^ x + EE 7 ^ x d k ^ 

(12) 

pm,, > N m v ,t (13) 

2.1.3. Facility allocation module 

The costs related to CSP (Z p ) and biorefinery facilities (Z r ) are 
composed of annual operating and capital related costs (Eq. (14) 
and (15)). Annual operating costs include the costs for utilities, 
maintenance, labor, supervision, insurance, laboratory charges, 
and waste treatment. In this study, it is assumed that the facilities 
with different capacities incur the same unit operating costs {p p 
and Pr). Therefore, annual operating costs are linearly dependent 
on the amount of biomass processed. Facility capital related costs 
are estimated by considering fixed (if p{ and ) and variable (i)^ 
and if Tp ) capital costs. Because of the economies of scale, fixed 
and variable facility capital costs vary at different ranges of 
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capacity. A piece-wise linear approximation approach was previ¬ 
ously developed to estimate the associated capital costs at differ¬ 
ent capacity scales (Lin et al., 2013) and applied to this study. 
The key decisions for the facility module include facility capacities 
(Cp S and Cj. s ) and locations (0* and 0() on the strategic level and 
facility operating capacity (f* f and P'f) and inventory monitoring 
(/*'*) on the tactical level. 

z *=EE(° p s >< ** + x «*,) + EE^ x (i4) 

4 = EE (°* x % + x Of.) + ^ x n r (15) 

The amount of biomass received at a CSP facility at each time 
period (A^ f ) can be estimated based on the daily delivery schedules 
from all farms (D IJ ' kt ) considering biomass loss during transporta¬ 
tion and handling (a>„) (Eq. (16)). The integrated model assumes 
that all the biomass feedstock received at a facility would be stored 
in a biomass input storage site first and then be handled for pre¬ 
processing. Given the differences between the amounts of biomass 
received and processed (P£ f ) at the each time period, the amount of 
unprocessed biomass feedstock inventory (Q p f ) could be moni¬ 
tored at each time period considering the unprocessed biomass 
storage loss (oj Pi ) (Eq. (17)). The designed facility capacity (C^) 
should be greater than or equal to the amount of biomass pro¬ 
cessed at any time period (Eq. (18)). A piece-wise linear approxi¬ 
mation is applied to estimate the economies of scale for CSP 
facilities (Eq. (19)-(21)), and similarly for biorefineries. The sum 
of the capacities at each capacity level (C p s ) equals the CSP facility 
capacity for that county (Eq. (19)). The facility capacity at each 
capacity level is constrained by the capacity range of each level 
(2| and /. p ) and its capacity level selection (0£ ,s ) (Eq. (20)). If there 
exists a facility in county k, the binary variable Op will be set equal 
to one, and exactly one binary variable O ks will be allowed to equal 
one (Eq. (21 )). Alternatively, if no facility exists in county k, O k will 
be set equal to zero, and all binary variables (Op’ s ) must be equal to 
zero (Eq. (21 )). Preprocessed biomass inventory at each CSP facility 
by the end of each time period (l p l ) can be estimated by the accu¬ 
mulated differences between the amounts of biomass processed 
and biomass delivered to all biorefineries (D IJ fct ) taking into con- 


sideration the processed biomass storage loss (a> P2 ) (Eq. (22)). 

EE D ' mt 

X (1 - to„) = A k / 

(16) 

<#*-(<£*■ 

1 + (Ap 1 - if) x At) x (1 - cop, x A t ) 

(17) 

pk,t < £lt 


(18) 



(19) 

4 X 0 k / 4 c k / 4 4 X oM 

(20) 

E o k / = 0 p 


(21) 

i k /=Ur 

+ (V/-ED k ’ , ’ t ] xAt] x(l-« P2 xAt) 

(22) 


This study assumes the biorefinery adopts just-in-time opera¬ 
tions, where the amount of biomass processed at a biorefinery 
facility at any time period (Pj f) is estimated based on the biomass 
delivery schedules from all CSP facilities and accounting for bio¬ 
mass loss by transportation and handling (Eq. (23)). The biorefin¬ 
ery processing throughput at any time period should not exceed 
its design capacity (Cj.) (Eq. (24)). For each biorefinery, the amount 


of ethanol produced (/i is the biomass-ethanol conversion rate) 
should be greater than or equal to the amount of ethanol trans¬ 
ported from its facility to all blending locations (Eq. (25)). 

E 0 *’ 4 x (1 - oi v ) = P 1 / (23) 

P'f <: C' (24) 

E^x^EED'^ ( 25 ) 


2.1.4. Ethanol distribution module 

Ethanol distribution costs (Z d ) are estimated based on the unit 
ethanol transportation cost (jua\ ethanol transportation patterns, 
and the distance between biorefineries and blending stations as 
well as between blending stations and ethanol consumption sites 
(Eq. (26)). The model assumes ethanol producers would outsource 
ethanol distribution to third-party logistics companies based on a 
unit transportation cost, in terms of $ gar 1 km -1 . Ethanol distribu¬ 
tion patterns between ethanol plants and blending stations (D lb ‘) 
and between blending and gas stations ( D bet ) are the key decision 
variables, whereas unit ethanol transportation cost (/t d ), transpor¬ 
tation distances ( e'' b and 0 bc ), and county-level ethanol consump¬ 
tion demand (£') are key input parameters. 

For each ethanol blending station, the total amount of ethanol 
received from all biorefineries should be greater than or equal to 
the amount of ethanol distributed to all possible ethanol consump¬ 
tion locations (Eq. (27)). For each ethanol consumption site, the to¬ 
tal amount of ethanol sourced from all possible blending stations 
should meet its ethanol consumption demand (E 6 ) (Eq. (28)). 

4 = (EEE D ' At x o ,b +EEE Di, ’ et x ^ x Md (26) 
EE d ' ,m > EE D ' ,,et (27) 

EE D ' , ’ e ’ t = F (28) 


2.2. Workflow of the integrated model 

Based on our preliminary tests, the integrated strategic and 
tactical optimization model requires significant computational 
resources and could become challenging for large-scale systems 
optimization. Therefore, a two-step workflow was proposed in or¬ 
der to develop an efficient approach to optimize a large-scale bio¬ 
fuels production system (Fig. 2). The strategy behind this approach 
is first to select a sub-optimal supply chain region using a strategic 
planning model, and then to optimize both strategic and tactical 
planning decisions simultaneously within the sub-optimal region 
using the integrated optimization model. 

Given an analysis region defined by users, county-level biomass 
production costs and biomass availability will be estimated first 
via the BioFeed model (Shastri et al., 2010). Biomass yield, farm 
size, and equipment data are the key input parameters for farm 
production optimization. The estimated biomass production costs 
and availability, together with system supply and demand require¬ 
ment, will be passed as inputs to the previously developed strate¬ 
gic model, the BioScope model (Lin et al„ 2013). The BioScope 
model will be executed to select the optimal biomass supply chain 
region to meet the defined biomass demand. Since the BioScope 
model was developed to optimize long-term planning decisions, 
only general estimations of biomass supply, processing, and 


T. Lin et al. /Bioresource Technology 156 (2014) 256-266 


261 



Fig. 2. A two-step workflow of strategic and tactical supply chain optimization for large-scale biofuels production. The first step is to select a sub-optimal biomass supply area 
to meet the demand, and the second step is to optimize both strategic and tactical decisions within that sub-optimal area. Rectangles are input parameters, rounded 
rectangles are models, and hexagons are the results from the integrated model. 


storage performance are considered in the model. To ensure that 
the selected region would meet the targeted biomass demand in 
the detailed analysis, a conservative estimate of biomass provision 
performance, such as biomass loss rate during biomass harvesting, 
transportation, and storage, will be applied in the BioScope model. 
As a result, the selected biomass supply area from the BioScope 
model should require larger areas than what the true optimal solu¬ 
tion would select. The selected sub-optimal area will be passed as 
the possible biomass supply counties for the integrated strategic 
and tactical optimization in the second step. The strategic and tac¬ 
tical decisions will then be optimized simultaneously within that 
sub-optimal area, based on the constraint equations listed in the 
previous section. 

Through this two-step workflow, the selected potential biomass 
supply chain region after the first step could be much smaller than 
the initial user defined region. The reduced number of possible 
supply counties and farms would save significant computational 
effort in solving the integrated strategic and tactical decisions at 
farm level. The model is capable of providing logistics decisions 


from a daily basis to a weekly basis in accordance with the user’s 
requirement. 

2.3. Case study for model applications 
2.3.1. Base case scenario 

To illustrate the use of the integrated biomass supply chain 
optimization model, we chose a Miscanthus-ethanol supply chain 
in Illinois for a case study. In the baseline case, 5% of cropland at 
each of 102 counties in Illinois was assumed to be allocated for 
Miscanthus production. Given the number of farms and their crop¬ 
land size distributions (USDA, 2010), 5% of farms in each size cat¬ 
egory in each county would be considered as potential biomass 
supply sites. The number of possible farms in each county of Illi¬ 
nois ranged from 4 to 82. The largest city in each county of Illinois 
was considered a candidate location for building CSP and biorefin¬ 
ery facilities, where either stage had 102 candidate locations. The 
base case study was designed to produce sufficient ethanol to 
replace 1% of gasoline consumption in the eight counties of the 
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Greater Chicago area. The county-level gasoline consumption data 
are based on county-level populations (US Census Bureau. 2012) 
and per capita gasoline consumption (US EIA, 2011). The total eth¬ 
anol requirement was approximately 32 million gallons per year 
for 1% gasoline replacement in the Greater Chicago area (Supple¬ 
mentary Fig. S.l ). Assuming a biomass-ethanol conversion rate of 
80 gal Mg 1 , 400,000 Mg of biomass was required for conversion 
annually. Among these eight counties, four counties had existing 
blending facilities (Supplementary Fig. S.l). 

Miscanthus production in Illinois reaches its peak yield in win¬ 
ter and peak yield data in each county were based on the MI5CAN- 
MOD model (Jain et al., 2010). Miscanthus was assumed to be 
harvested from January to April in this study. During the harvest¬ 
ing season, the Miscanthus dry matter yield was reduced by 
approximately 0.07 Mg ha -1 d 1 from the peak yield, while the 
moisture content was reduced by 0.1% d 1 from 25% to a minimum 
of 15% (Shastri et al„ 2010). The probability of a working day dur¬ 
ing the harvesting season varied with time as well (Shastri et al., 
2012). Given the existing interstate and state highway road net¬ 
work (US Census Bureau, 2011), ArcGIS was used to calculate the 
shortest distance between any potential biomass supply, CSP, 
and biorefinery. Due to the lack of specific farm location informa¬ 
tion, the transportation distance within each county from a farm 
to the largest city was approximated as a stochastic number from 
zero to the radius of a circle with the same area as the county. 

Farm management includes the following unit operations: har¬ 
vesting, baling, in-field transportation, and stacking. Biomass was 
stacked in bales at the farm-gate for delivery. Trucks were used 
to serve both stages of biomass transportation, from farms to CSP 
facilities and from CSP facilities to biorefineries. The capital costs 
and operating data of farm equipment and vehicles were adopted 
from Shastri et al. (2010). CSP facilities were designed to reduce 
biomass particle size through tub-grinding. The ground biomass 
was further tapped to increase the density to 200 kg m 3 , higher 
than 150 kg nrr 3 of baled biomass (Sokhansanj et al., 2009). 

Due to the economy of scales, the facility capital investment 
costs for biorefineries and CSP facilities was estimated using a scal¬ 
ing factor of 0.7 (Lin et al., 2013). The baseline case costs were $422 
million for a biorefinery at a capacity of 2144 Mg d 1 (Humbird 
et al., 2011) and $19 million for a CSP facility at a capacity of 
2016 Mg d 1 (Hess et al., 2007). An annualized cost factor of 
13.7% was used to estimate the annualized capital related costs 
(Lin et al., 2013). The capacity level and its associated capital and 
operating costs for both CSP and biorefinery facilities were based 
on Lin et al. (2013). All the costs used in this study have been con¬ 
verted to year 2007 followed the procedure previously employed 
in the strategic modeling analysis (Lin et al. 2013). 


2.3.2. Scenario analysis of Miscanthus yield change 

As would be expected with crop growth, Miscanthus yield is 
highly dependent on moisture stress, which varies from season 
to season. A sensitivity analysis was conducted to illustrate impact 
of changes in Miscanthus peak yield on the optimal biomass supply 
chain configuration. Especially, the model could quantify how the 
change of biomass yield would affect farmers’ tactical operating 
decisions. The county-level Miscanthus peak yield was considered 
at 100%, 75%, and 50% of the data from the MISCANMOD model 
(Jain et al., 2010) for each county in Illinois. 5% of cropland was 
allocated for producing Miscanthus and the ethanol demand was 
targeted for 1% of gasoline replacement in the Greater Chicago Area 
in all three scenarios. 

2.3.3. Scenario analysis of demand change 

A sensitivity analysis was conducted to illustrate the impact of 
the change of biomass demand on the optimal biomass supply 
chain configuration. Three levels of ethanol demand were consid¬ 
ered at 1%, 2%, and 4% of gasoline consumption replacement in 
the Greater Chicago Area, respectively. The required biomass de¬ 
mand increased from 400,000 to 1,600,000 Mg per year at 4% gas¬ 
oline replacement rate. 5% of cropland was allocated for producing 
Miscanthus in all three scenarios. 

3. Results and discussion 

3.1. Base case scenario 

Two counties (Iroquois and Kankakee counties) were selected to 
support 32 million gallons of ethanol production. The model sug¬ 
gested building one centralized storage and preprocessing (C5P) 
facility and one biorefinery in Kankakee County. Although biomass 
yield is relatively high in southern Illinois counties, the selection of 
these two counties, which are located near the Chicago ethanol 
consumption area, is the result of balancing the tradeoff between 
biomass procurement costs and ethanol distribution costs. The 
optimal Miscanthus-ethanol production costs were $218.5 Mg -1 
of biomass, or $0.72 L 1 ($2.73 gal -1 ) of ethanol. Biorefinery 
related costs were the largest cost component, accounting for 62% 
of the system costs, followed by biomass procurement, C5P 
related costs, biomass transportation costs, and ethanol distribu¬ 
tion costs (Table 2). Because biorefinery related costs have a signif¬ 
icant impact on the total costs, it is worthwhile to build a 
large-scale facility to gain the economies of scale. The average 
biomass procurement costs were $52.8 Mg -1 , which is in the previ¬ 
ously observed range of $41 to 54 Mg -1 (Khanna et al„ 2008). The 
unit biomass production cost (excluding the land opportunity 
costs) was $23.1 Mg~\ where baling is the most expensive farm 


Table 2 

The comparison of unit Miscanthus-ethanol production costs at different Miscanthus yield rates. The base case Miscanthus yield is based on Jain et al. (2010), where the county- 
level peak yield ranges from 38 to 48 Mg ha 1 in Illinois. All three scenarios assume 5% cropland is allocated to grow Miscanthus to support 1% of gasoline replacement in the 
Greater Chicago area, an annual demand of 32 million gallons of ethanol. All the numbers have been converted into $ per Mg of dry biomass, assuming a biomass-ethanol 
conversion rate at 80 gallons per Mg of dry biomass. 


Biomass procurement costs (A) 

-Harvesting costs 
-Baling costs 

-Infield transportation costs 
-Handling costs 
-Land opportunity costs 
Transportation costs (B) 

CSP facility related costs (C) 

Biorefinery related costs (D) 

Ethanol distribution costs (E) 

Unit Miscanthus-ethanol production costs (A +1 


Base case scenario Miscanthus ; 


52.8 60.6 

3.1 3.3 

15.5 15.2 

3 3.4 

2.3 2.6 

28.9 36.1 

11.8 10.5 

14.7 15.8 

135.8 135.8 

3.4 7.7 

C + D + E) 218.5 230.4 


: 75% of base case Miscanthus yield at 50% of base case 


4 

14.3 

3.1 

2.2 
40.8 

14.7 

135.8 

16.5 

242.5 
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County level data regarding the number of farms selected to produce biomass, their area, biomass average deliverable yield, capture rate, and the amount of biomass supplied. 

Biomass supply Number of farms selected in Amount of selected cropland Average deliverable yield Biomass capture rate (as of Supply amount 

county each county area (ha) (Mgha -1 ) peak yield) (Mg) 

Iroquois 23 8779 28.3 80% 248,600 

Kankakee 17 6584 26.9 80% 176,860 


operation, followed by harvesting, in-field transportation, and han¬ 
dling. The cost is lower than a cost of $31.2 Mg 1 from the baseline 
tactical model (Shastri et al, 2010) and $40 Mg -1 in Khanna et al, 
(2008), which possibly results from the large farm selection and 
high deliverable yield of Miscanthus in the current study. Biomass 
land opportunity costs, $28.9 Mg~\ accounted for more than half 
of the procurement costs. 

Within the two selected counties, 40 out of 126 possible farms 
were selected to produce 425,000 Mg of Miscanthus annually 
(Table 3). The total biomass losses during storage, preprocessing, 
and transportation steps account for 25,000 Mg, or 6%. All the se¬ 
lected farms are large farms due to their relatively low biomass 
production costs, where the size of farms ranged from 180 to 
1610 ha, or 445 to 3977 acres. The deliverable biomass yield ran¬ 
ged from 29.6 to 31 Mg ha~\ which represents 80% of peak yield. 
The average deliverable yields of these two counties selected are 
higher than a state-level average yield of 22.3 Mg ha -1 in Illinois 
(Khanna et al., 2008). The biomass loss on farms mainly results 
from three parts: standing biomass dry matter loss in the field, 
equipment operating loss, and storage loss. The number of pieces 
of farm equipment required for each farm varied significantly 
given the differences in their cropland areas. The largest farm se¬ 
lected in this study could provide 34,752 Mg of biomass annually. 
This farm required 5 harvesters, 14 balers, 8 in-field truck trailers, 
and 4 loaders. The large number of balers was mainly due to the 
low throughput rates assumed in the model. 

The proposed delivery schedule recommended delivering more 
biomass from farms to the C$P facility in the biomass harvesting 
season (first 120 days), to reduce biomass losses on farm fields. 
The C5P facility operated at a constant rate of 1136 Mg d \ The 
C5P biomass input rate was much larger than its preprocessing rate 
in the first 120 days, which resulted in an accumulation of biomass 
feedstock inventory (Fig. 3). The biomass feedstock inventory 


peaked at the end of harvesting season with the amount of 
248,285 Mg, which requires 43 hectares, or 106 acres, of land for 
biomass feedstock storage. The accumulated baled biomass would 
be further preprocessed during the last 240 days of the year to 
maintain the design processing rate. The total annual operating 
capacity of the C5P facility was 408,915 Mg, which was 3.7% lower 
than the total supply of biomass from all farms at 425,000 Mgs. The 
3.7% biomass loss is a result of preprocessing operation loss and 
biomass feedstock storage and handling loss at C5P facilities. 

The delivery rate of preprocessed biomass (ground biomass 
with compaction) is constant throughout the year. The system 
achieved just-in-time biomass delivery for a biorefinery operating 
at a daily operating capacity of 1102 Mg d '. With the manage¬ 
ment of logistics planning, the system required 54 vehicles to fulfill 
the delivery schedules. Among them, 49 vehicles were assigned for 
delivery between farms and C5P facilities, as a result of distributed 
locations of biomass supply farms. 


3.2. Miscanthus yield change impact 

The optimal biofuel production costs increased to $242.5 Mg -1 
from $218.5 Mg~\ or an 11% increase, if the Miscanthus peak yield 
would be 50% of that in the base case. The key factor is the increase 
of biomass procurement and ethanol distribution costs (Table 2). 
The significant increase of the biomass procurement costs is largely 
due to the increase of the unit land opportunity costs of growing 
Miscanthus. In order to convince farmers to grow Miscanthus, 
the model assumes that farmers would gain at least the same profit 
per unit area as producing com and soybean on the same land 
(Khanna et al., 2008). The decreased biomass yield, therefore, 
would result in the increased land opportunity costs per unit bio¬ 
mass procured. 



Fig. 3. The optimiz 


CSP operating information. 
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County level data regarding the number of counties and farms selected to produce biomass, their area, and biomass average deliverable yield at different levels of biomass yield 
and ethanol demand. 


Scenario analysis 


Ethanol demand as of gasoline replacement rate 


Number of selected Number of farms Selected cropland Average deliverable 

counties selected area (ha) yield (Mg ha~’) 


Base case yield 1% 

75% of base case yield 1% 

50% of base case yield 1% 

Base case yield 2% 

Base case yield 4% 


8 


40 15,363 27.7 

37 18,425 23.1 

58 23,453 18.1 

83 27,413 28.9 

203 57,450 27.6 


Table 5 

The comparison of unit Miscanthus-ethanol production costs at different ethanol demand levels. All three scenarios assume 5% cropland is allocated to grow Miscanthus. Ethanol 
demand of the base case scenario is to replace 1% of gasoline demand in the Greater Chicago area that requires an annual demand of 32 million gallons of ethanol. All the numbers 
have been converted into $ per Mg of dry biomass, assuming a biomass-ethanol conversion rate at 80 gallons per Mg of dry biomass. 

Ethanol demand 


Biomass procurement costs (A) 

-Harvesting costs 

-Infield transportation costs 
-Handling costs 
-Land opportunity costs 
Transportation costs (B) 

CSP facility related costs (C) 

Biorefinery related costs (D) 

Ethanol distribution costs (E) 

Unit Miscanthus-ethanol production costs (A + B + C + D + E) 


1% Gasoline replacement 2% Gasoline replacement 4% Gasoline replacement 


52.8 

3.1 

15.5 

3 

2.3 

28.9 


135.8 

3.4 

218.5 


51.7 
3.1 

15.9 
3.1 

2.9 
26.7 

11.7 
14.2 
123.5 

7.7 
208.8 


52.8 

3.7 
15.3 
3 

2.7 
28.1 

15.2 

108.8 
5.3 
195 


In all three scenarios, one centralized biorefinery was suggested 
to gain the economies of scale. With the decrease of biomass yield, 
the selected biomass supply counties were moved to southern Illi¬ 
nois, where relatively high biomass yields were still available (GIS- 
enabled result visualization is provided in Supplementary Fig. S.2). 
The distance of moving ethanol has increased from less than 80 km 
to more than 500 km at the low biomass yield case. As the ethanol 
facility moves away from Chicago area, the number of selected 
blending stations increased from two to three for better ethanol 
distribution. The results showed that the increase of biomass pro¬ 
curement costs outweigh the savings of ethanol distribution costs 
if biomass supply remained near Chicago area. At low biomass 
yield rate, the system would choose a high biomass yield area as 
its supply, located in Southern Illinois. 

The required number of biomass supply counties increased 
from two to five, when biomass yield rate reduced to half of the 
base case (Table 4). The number of farms increased to 58 from 
40 in the base case, with larger average farm size in the Southern 
Illinois counties. Although the projected biomass yield at each 
county was reduced to half of the base case, the required produc¬ 
tion area increased to 23,453 from 15,363 ha, representing only a 
50% increase of the base case, due to the relative high biomass 
yield in Southern Illinois. 

3.3. Miscanthus demand change impact 

The total demand of biomass for conversion increased from 
400,000 to 1,600,000 Mg annually at a 4% gasoline replacement 
rate. The optimal biofuel production costs were reduced with high¬ 
er annual biomass demand, from $218.5 Mg 1 of base case to 
$195 Mg -1 at an annual demand of 1.6 million Mg (Table 5). In 
all three scenarios, the system suggested building a centralized 
biorefinery to gain the economies of scale, while the number of 
biomass supply counties and CSP facilities increased with the 
demand increase. Because of the economies of scale, the biorefin¬ 
ery related costs could be reduced by about $27 Mg -1 if biorefinery 


capacity was four times that of the base case. The savings in bior¬ 
efinery related costs exceeded the increases in biomass procure¬ 
ment and transportation costs as a result of a larger biomass 

With higher demand, the system suggested a distributed bio¬ 
mass supply chain configuration with more C5P facilities (GI5- 
enabled result visualization is provided in $upplementary 
Fig. $.3). The system required 203 farms from eight counties to 
support 4% gasoline replacement, whereas 40 and 83 farms were 
required for 1% and 2% gasoline replacement, respectively (Table 4). 
The total required biomass production area increased from 15,363 
to 57,450 ha when the ethanol demand increased to replace 4% of 
gasoline consumption. The fleet size also increased from 54 vehi¬ 
cles for the base case to 156 vehicles for the largest demand case. 

3.4. Model limitations and future works 

The major limitations of the present integrated model are data 
availability and uncertainty. Detailed farm-level data, including 
farm locations and sizes, and the options of farming equipment 
and associated cost and operating performances, are needed for a 
thorough farm management analysis. Furthermore, biomass pro¬ 
duction and probability of working days change annually given 
weather uncertainties. $tochastic optimization approach could be 
incorporated into the present model to better understand the 
uncertainty issues related to biomass production. The other major 
limitation is that the present model assumes biorefinery compa¬ 
nies control the whole supply chain, providing biomass delivery, 
storage, and preprocessing services. Future research will be con¬ 
ducted to incorporate agent-based modeling approach to under¬ 
stand the behaviors and interactions among different parties in 
the system. 

The integrated strategic and tactical planning model developed 
in this study could be applied to facilitate biomass provision con¬ 
sidering different feedstocks such as other types of energy crops 
and agricultural residues in different regions. The model would 
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also be used to evaluate different preprocessing and conversion jk,t 
technologies (Dwivedi et al., 2009; ElMekawy et al., 2013) by incor- p 

porating detailed processing information. For the long term, it is 
expected that the supply chain optimization model could be inte- jjijM 
grated with environmental and life cycle analysis (Singh et al., 

2010) to evaluate the environmental and economic performance 
ranging from regional specific feedstock production to detailed 
biomass conversion processing. 


4. Conclusion 

An integrated strategic and tactical supply chain optimization 
model was developed to minimize annual biomass-biofuel pro- 
duction costs. By considering spatial and temporal constraints, 
the model can identify how long term decisions, such as facility F™" 
locations and capacities, interact with short-term decisions, such qKs 

as weekly biomass production and delivery schedules. Several case p 

studies of Miscanthus-ethanol production in Illinois were pre- A k,t 

sented to illustrate the usage of the model. The results of scenario p 

analyses showed that the changes of biomass yield and ethanol 
demand caused a significant impact on costs and supply chain pk,t 
configurations. p 
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Appendix A. Nomenclature 


Set Set name and element labels 

/ Biomass supply county, where i e I 

J Farm, where j e] 

K CSP facility location, where k e I< 

L Biorefinery location, where l e L 

B Blending location, where b e B 

E Ethanol consumption county, where e e E 

T Time period, where t eT. t*' is an element 

representing the end of the harvesting season 
and f is an element representing the end of 
the simulation cycle 

M Type of equipment, where m e M. m f is a 

subset of M representing farming equipment 
and m v is a subset of representing road 
vehicles. 

S Capacity scale for a facility, where s e S 

Decision Description 

variables 

pi-j.tnf The number of each type of equipment 

f required for farm j in county i 

p'J.m/.f The operating hours of each type of 

fs equipment required for farm j in county i at 

time period t 

jfit The amount of biomass delivered to the gate 

of farm j in county i at time period t (Mg d 1 ) 
H'it The cumulative harvested area for farm j in 

county i until the end of time period t (km 2 ) 
fit The biomass inventory for farm j in county i 

at the end of time period t 


C* 


D l,bt 


°P 

°r 

Of 


Inputs 


The preprocessed biomass inventory for CSP 
facility located in county k at the end of time 
period t 

The amount of biomass delivered from farm j 
in county i to CSP facility k at time period t 
(Mgd- 1 ) 

The number of trips required from farm j in 
county i to CSP at county k by vehicle m v at 
time period t 

The number of trips from CSP in county k to 
biorefinery in county / by vehicle m v at time 
period t 

The number of vehicle type m v required at 
time period t 

The required fleet size of each type of vehicle 
The capacity of CSP facility in county k at the 
capacity scale s (Mgd ') 

The amount of biomass received by CSP 
facility in county k during time period t 
(Mgd- 1 ) 

The amount of biomass preprocessed by CSP 
facility in county k during time period t 
(Mgd- 1 ) 

The unprocessed biomass inventory for CSP 
facility located in county k at the end of time 
period t 

The amount of biomass delivered from CSP 
facility in county k to biorefinery in county / 
at time period t (Mg d 1 ) 

The amount of biomass processed by 
biorefinery in county / at time period t 
(Mgd- 1 ) 

The capacity of biorefinery in county / 
(Mgd- 1 ) 

The amount of ethanol delivered from 
biorefinery in county / to blending station in 
county b at time period t (Mg d 1 ) 

The amount of ethanol delivered from 
blending station in county b to ethanol 
consumption site in county e at time period t 
(Mgd- 1 ) 

Binary variable: Whether farm j in county i is 
selected for biomass production 
Binary variable: Whether a CSP facility is 
located in county k 

Binary variable: Whether a biorefinery is 
located in county / 

Binary variable: Whether a CSP facility is 
located in county k at the capacity scale s 
Binary variable: Whether a biorefinery is 
located in county 1 at the capacity scale s 

Description 

The unit operating cost of each type of 
equipment ($ h ’) 

The annual capital related cost of each type of 
equipment 

The size of each farm j in county i (km 2 ) 

The unit cropland opportunity cost in county 
i ($ km" 2 ) 

Operating hours per day 
The travelling time from farm j in county i to 
CSP facility in county k during time period t 
using vehicle type m v (h) 


(continued 


page) 



?? ? | | 8,8 .* & ^ ^ % * 


d k ,i,m„ The travelling time CSP facility in county k 

during time period t to biorefinery in county / 
using vehicle type m v (h) 

The weight limit of each type of vehicle (Mg) 
The volume limit of each type of vehicle (m 3 ) 
The annual capital related cost of each type of 
vehicle 

The unit operating cost of each type of 
vehicle ($ h 1 ) 

The unit biomass handling cost at CSP 
facilities ($ Mg 1 ) 

Ethanol conversion rate (gal Mg -1 ) 

The distance between biorefinery location in 
county / to ethanol blending station in county 
b (km) 

The distance between blending station in 
county b and ethanol consumption in county 
e (km) 

The density of delivered biomass from farm j 
in county i (Mg nT 3 ) 

The fixed capital cost to build a CSP facility at 
the capacity scales 

The unit variable capital cost to build a CSP 
facility at the capacity scale s ($ Mg -1 ) 

The fixed capital cost to build a biorefinery at 
the capacity scale s 

The unit variable capital cost to build a 
biorefinery at the capacity scale s ($ Mg -1 ) 
The biomass loss rate occurred at farms 
The biomass loss rate occurred during 
transportation and handling 
The unprocessed biomass loss rate during 
storage at CSP facilities 
The processed biomass loss rate during 
storage at CSP facilities 
Unit operating cost of ethanol production 
(SMg- 1 ) 

Unit operating cost of preprocessing ($ Mg" 1 ) 
The unit ethanol distribution cost 
($ gal" 1 km- 1 ) 

[i hl The unit biomass handling cost at 


Appendix B. Supplementary material 

Supplementary data associated with this article can be found, in 
the online version, at http://dx.doi.org/10.1016/j-biortech.2013. 
12.121. 















